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O ' Abstract. White-light observations of the solar corona show that there are two character- 

H " 

^ ' istic types of Coronal Mass Ejections (CMEs) in terms of speed-height profiles: so-called 

fast CMEs that attain high sp eeds low in the corona and slow CMEs that gradually accel- 



erate from low initial speeds. 



Low & Zhang 



2002h have recently proposed that fast and 



slow CMEs result from initial states with magnetic configurations characterized by nor- 
mal prominences (NPs) and inverse prominences (IPs), respectively. To test their theory, 
we employed a two-dimensional, time-dependent, resistive magnetohydrodynamic code to 
simulate the expulsion of CMEs in these two different prominence environments. Our nu- 
merical simulations demonstrate that (i) a CME-like expulsion is more readily produced 
in an NP than in an IP environment, and, (ii) a CME originating from an NP environment 
tends to have a higher speed early in the event than one originating from an IP environ- 
ment. Magnetic reconnection plays distinct roles in the two different field topologies of 
these two environments to produce their c haracteristic CME spe ed-height profiles. Our nu- 



merical simulations support the proposal of lLow & Zhans 



2002» although the reconnection 



development for the NP associated CME is different from the one sketched in their theory. 
Observational implications of our simulations are discussed. 
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1. Introduction 



The accumulated observations of Coronal Mass Ejections (CMEs) obtained over more than 
3 solar cycles by the coronagraphs onboard the Skylab, Solar Maximum Mission, and Solar 
and Heliospheric Observatory (SOHO) have shown two characteristic types of CME speed- 
height profiles: (I) so-called fast CMEs attaining high speeds above the CME median speed 



(II) slow CMEs gradually accelerating from low initial speeds (< 400 km s ^) i 


Goslint 


' et al, 


1976 


MacOueen & Fishei 


1983; 


Drvei 


1996; 


Sheelev et al. 


1999; 


St. Cvr et al. 


I999I 


2000; 



Moon et al J2OO3. t has also been recognized that Type I CMEs tend to originate from active re- 
gions and are frequently accompanied by flares, while Type II events usually originate away from 
active regions and are accompanied by eruptive prominences. Outstanding examples of these two 
characteristic types can be found in a given sample of CMEs, but it should be pointed out that 
there is a continuous spectrum of speed-height profiles between these two types as the opposite 



extremes. Another way to distinguish these two types of events was given bv lAndrews & Howar( 



] 



1 20011) . They showed that CMEs observed by the LASCO coronagraphs onboard SOHO may fall 
into two categories: Type C characterized with a constant speed and Type A characterized with a 
significant acceleration, which respectively correspond to Types I and II. 

The underlying physics responsible for the dual character of CME speed-height profiles re- 
mains among the outstanding questions in CME research. Recently, Low & Zhang (2002, here- 
after referred to as LZ02) suggested a theory to explain this phenomenon in terms of the different 
hydromagnetic environmen ts in which CMEs occur. CMEs are correlated at about the 75% leve l 



with prominence eruptions ( iMunro et alJll97S ; 



Webb & Hundhausenll987 



I 



St. CvretalJl999h 



There are tw o magnetic types of 



quies cent prominences referred to as the Normal and Inverse 



Prominences dTandberg-Hanssenl 1 9951) . hereafter called NPs and IPs respectively. Low & Zhang 
pointed out that these two types of prominences represent different magnetic field topologies 
which have distinct consequences for the interplay between magnetic reconnection and CME ex- 
pulsion dynamics. They gave qualitative sketches of this hydromagnetic interplay which suggest 
that fast and slow CMEs are naturally associated with NPs and IPs, respectively. 

In the case of an NP, the surrounding field has a topology such that the CME is expelled with 
a current sheet to be dissipated by magnetic reconnection ahead of the erupt ing prominence. 



This re connection is a break-out effect similar to the one originally proposed by 



Antiochos et al. 



1 19991) for a multipol ar magnetic field. In the LZ0 2 proposal, the global field is bipolar with 



a magnetic flux rope dChenlll989 ; 



Chen et al 



19971) . The rise of the flux rope and prominence 



drives break-out reconnection ahead. Reconnection produces a slingshot effect which, in turn, 
drives the CME and prominence. This runaway situation naturally produces a fast CME with 
an impulsive acceleration and flare heating early in the event. It is worth mentioning that such 
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a timing correlation between a fast CME and flare does not necessarily mean that the flare is 
driving the CME. 

In contrast, an IP is associated with a magnetic-flux rope topology which produces a current 
sheet trailing behind the CME and erupting prominence. This current sheet is not driven directly 
by the CME but forms passively from the left-behind magnetic field. Magnetic reconnection is 
then not compelled to occur early and there is no impulsive CME acceleration or flare heating 
by reconnection early in the event. A CME produced in this hydromagnetic environment may be 
expected to be gradually accelerated from a low initial speed. 



St udies of flare morphol ogy and timing associated with fast and slow CMEs dZhang et al. 



2002; 



Zhang & Golut 



2003^ have shown that observations of this kind are consistent with the 
proposed theory of LZ02. Ultimately, observations of more incisive kind will be needed to verify 
or reject the LZ02 theory. On the other hand, the proposed theory of LZ02 is based on intuitive 
sketches of the relevant hydromagnetic processes. It is therefore important to perform numer- 
ical magnetohydrodynamic (MHD) simulations to directly investigate if the two distinct CME 
acceleration processes of LZ02 can be produced. 

Individual si mula tion works have b een reported in the literature on the acceleration of CMEs. 



Wu et al 



1 19951) and 



Wu et al 



( 1997h . respectively using a 2-D and 2.5-D ideal MHD model, 
studied CME expulsions involving a magnetic flux rope that can be interpreted to be those repre- 
senting an IP. In their models, there was no magnetic reconnection taking place. A si milar study. 



but inv olving reconnection in a configuration resembling an NP, was performed by 



Guo et al. 



1 1996h using a 2-D resistive code. However, the topological differences between the two types 
of configuration had not been recognized in connection with the two-class CMEs by these au- 
thors, nor by others until th e recent work by LZ02 . On the other ha nd, different numerical mod- 



els adopted in these works jWu et al 



1995 



1997 



Guo et alJI 19961) make it difficult to directly 



compare the CME acceleration between the two configurations from their results. Therefore a 
systematic, comparative simulation study of CME expulsions involving magnetic reconnection 
in both IP and NP configurations will be instructive in checking on the ideas of LZ02 and in an 
exploration of the interplay between magnetic reconnection and CME expulsion dynamics. Such 
a study is reported in this paper, and, as we shall see, new effects not considered by LZ02, with 
regard to CMEs in an NP environment, are among the results we will report. § 2 describes the 
numerical model. Simulation results are presented in § 3, followed by discussions on operating 
forces, magnetic topologies, and magnetic reconnection in § 4. Finally, we conclude this paper 
in § 5 with remarks relevant to the LZ02 theory and observational implications. 



2. Description of the Simulation IVIodel 

The n umerical model used for this study, based on those of 



Wu et al. 



1 19951) and 



Guo et al. 



1996I) . is composed of a set of two-dimensional, time-dependent, resistive, single-fluid MHD 



equations in the spherical coordinates (r, 6, (p), under an axisymmetry assumption (d/dcp - 0). 



W. Liu et al.: Magnetic topologies & two-class CMEs: a numerical MHD study 



The governing equations, including the conservation law s of mass, momen tum, and energy and 



Guo et al 



1 19961) which the interested 



the magnetic induction equation, are identical to those in 
reader is referred to for the mathematical formulation. We summarize the model as follows. 

(i) The computational domain is defined as IR^ < r < 7. 14/?, and -1.5° < 6 < 91.5°, where 
R, is the solar radius and - 0° (6 - 90°) the north pole (equator). We used an 81(r) x 63(6) 
grid, uniform in the meridional direction (AO = 1.5°) and non-uniform in the r-direction (Ar, = 
r, - = O.95[A0]r,_i) for a better resolution near the solar surface, with Ar, ranging from 
O.Q25Rs at the coronal base to 0.178/?.s at the outer boundary. The grid was staggered to prohibit 
sawtooth oscillations. 

(ii) The MH D equations were solved by the combined difference technique ( Guo et alii 991 



Wu et alJl995l) . The essence of this technique is to employ different numerical schemes to treat 



different equations according to their physical nature. Namely, the second-type upwind scheme 
is used for the continuity and energy equation, the Lax scheme for the momentum equation, and 
the Lax-Wendroff scheme for the magnetic induction equation. 

(iii) The boundary conditions are: (a) a symmetric boundary at the equator and pole, (b) 
linear extrapolation at the outer boundary, and (c) the method of projected normal characteristics 



I Wu & Wangll 987ft at the inner boundary. To guarantee the solenoidal condition, viz., V ■ B = 0, 
the reiterative divergence-cleaning method i Ramshawll983 ) was applied. 



(iv) One of the major modifications to the previous model was that we prescribed a uniform 
magnetic resistivity (rj = 8.75 x lO ^Q m) throughout the whole computational domain, rather than 



the anomalous resistivity used by 



Guo et al. 



19961) which favors reconnection only in regions 



with large electric current density. By doing this, various regions were treated equally, in view 
of the large difference in the current density distribution between the IP and NP configurations 
as we shall see in § 4. This modification is essential for a direct comparison between the two 
configurations. 

We too k three steps to si mulate the expulsion of a CME in the corona, following 



19951) and 



ep; 
ill 



IWu etal. 



Guo et alJ 1119961) . First we constructed an initial state of the corona with a quasi 



equilibrium helmet streamer by using the relaxation method dSteinolfson et al 



1982). The re- 



sulting magnetic field, plasma flow velocity, and electric current density are shown in Figure 1 
(see Fig. lb in Wu et al. 1995 for the corresponding plasma density distribution). The character- 
istic parameters are listed in Table 1 . 

In the next step, we emerged a flux rope with various energy contents and two types of mag- 
netic configuration from below the photosphere into the corona. Our 2-D model approximates the 
central cross-section of the 3-D flux rope which could be anchored at two ends on the photosphere 
in a reahstic geometry. The gas pressure and magnetic field of the flux rope (see Equations (4) 
and (5) in Guo et al. 1996), in an equilibrium state, were analytically specified in local cylindrical 
coordinates (r', ff ,z'), with Equation (5) being modified: 



(1) 
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where a is the radius of the flux rope. The "+" in this equation represents diff'erent senses of 
circulation (or polarity) of the flux rope's magnetic field: the poloidal field in the upper part of 
the rope with a "+" ("-") sign is of the same (opposite) circulation with respect to the external 
coronal field. Considering that the prominence material is presumably contain ed in the lower part 



of the flux rope, the "+" ("-") therefore corresponds to an IP (NP) topology ( iTandberg-Hanssen 



19951) . Equivalently, the azimuthal current in the emerging flux rope flows in the same direction 
for the "+" case as that in the vertical current sheet of the helmet streamer (Fig. lb). To implement 
the emergence process, we initially placed the flux rope below the photosphere with its center 
at r = Rs - a. We then slowly displaced the rope upward at a constant speed vvm va, see 
Table 2) by accordingly changing the physical variables at the inner boundary. It took 4 hours for 
the flux rope to entirely emerge into the corona and then the inner b oundary conditions were set 



Wu et al. 



1 19951) for a detailed mathematical 



back to their original form. The reader is referred to 
treatment of this flux rope emergence process. The last step was carried out by simply letting the 
system evolve by itself until f = 40 hrs. We present the simulation results in the next section. 



3. Simulation Results 

A total of 8 simulation cases, grouped into 4 pairs, were performed with various sizes of the 
emerging flux rope. The two cases of each pair have identical physical conditions except oppo- 
site magnetic polarities of the emerging flux rope, respectively corresponding to the NP and IP 
topology. The energy contents (i.e., the combined magnetic and thermal energy) inside the flux 
rope (assuming the third-dimensional thickness Az' = Q-lRs in the local cylindrical coordinates) 
are 1 .76 X 1 0^ 3 .95 X 1 0^ 7 .02 X 1 0^ and 1 . 1 X 1 0^^ ergs, respectively for the 4 pairs of cases 
(cf. ~ 10^^ ergs needed to expel a moderately large CME, see Hundhausen 1999; Forbes 2000). 
We tabulate the key parameters of these cases in Table 2 and describe detailed results as follows. 

Figure 2 shows the height-time and speed-time profiles of the "erupting flux rope"' center 
(defined as the O-type neutral point) for the studied cases. The solid (dotted) lines represent 
the cases of the NP (IP) environment. As we can see, the differences between the two types of 
environments are evident. In Cases la and 2a, the emerging flux rope has already destabilized 
the helmet streamer and launched a CME by the end of the simulation in the NP configuration. 
In contrast, the streamer and the flux rope are still in equilibrium in the IP configuration in Cases 
lb and 2b. As expected from LZ02's theory, these profiles reveal distinct characteristics between 
the two types of topologies: (i) an NP environment seems to be more in favor of producing 
CMEs than an IP environment (also see Zhang & Low 2003); (ii) a CME produced in an NP 
configuration tends to have a higher speed in its early life than a CME (if any) originating from 
an IP environment under otherwise identical conditions. We take note of that the average CME 

' For the NP configuration, the "erupting flux rope" refers to the new flux rope formed by reconnection; 
for the IP configuration, it refers to the originally emerging flux rope. This will be further explained in the 
following text. 
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eruption speed in Case 4a is about 2.3 times that in Case 4b (see Table 2), which agrees well 
with the typical fast-to-slow speed ratios of two-class CMEs derived from the Skylab and recent 
SOHO/LASCO coronagraphs (e.g. 775 km s"' versus 330 km s"' in Gosling et al. 1976, and 955 
km s ' versus 411 km s ' in St. Cyr et al. 1999). The other cases in this study show much higher 
ratios. 

Let us now focus our attention on two eruptive cases — 4a (NP configuration) and 4b (IP 
configuration). Figure 3 shows the magnetic field and the plasma velocity for these two cases. The 
corresponding azimuthal electric current density, J^, and plasma density enhancement, [p(r, 6, f)- 
p{r, 9, 0)]/p(r, 6, 0), are respectively displayed in Figures 4 and 5. 

3.1. Eruptive Case — Normal Prominence Configuration 

In Case 4a, we notice that upon the flux rope emergence, a curved current sheet is developed 
(Fig. 4a) between the leading edge of the emerging rope and the helmet dome (i.e., the closed 
field region of the helmet streamer) whose bipolar magnetic field is tied on the photosphere and 
in a direction opposite to the field in the upper part of the flux rope. Due to the finite resistivity, 
magnetic reconnection occurs in the current sheet and rapidly forms a new flux rope between the 
old (emerging) flux rope and the helmet dome (Fig. 3a). Once this new rope is formed, the current 
sheet splits into two halves (Fig. 4b) and two X-type neutral points appear on the flanks of the 
old rope (Fig. 3b). The magnetic flux in the new rope then grows as reconnection proceeds until 
aU the closed-field flux in the helmet dome is converted to the new rope's flux. In the meantime, 
the new flux rope expands and runs upward very fast while the old one follows and loses an equal 
amount of magnetic flux as the helmet-dome field above (Fig. 3b and 3c). Note that the old flux 
rope must contain sufficient poloidal flux to annihilate the closed field of the helmet dome ahead 
in order to set free the flux rope itself (also see LZ02). We identify the new flux rope as the main 
body of a CME in this numerical experiment, v/ and v in Table 2 respectively refer to the final 
and average speed of the new flux rope's center for the cases of the NP configuration. At ? = 25 
hrs, the erupting material has escaped from the computational domain with an average speed of 
V = 161.0 km s"^ and the system reaches a quasi-equilibrium state (Figs. 3d and 4d) similar to 
the initial state (Fig. 1). We further note that, during the eruption phase (i.e., t < 25 hrs), the 
maximal current is on the flanks of the old flux rope where the current sheets are located (e.g.. 
Figs. 4b and 4c) and the largest density enhancement occurs at the two lateral dips of the new 
flux rope, close to its center (the left column of Fig. 5). 

3.2. Eruptive Case — Inverse Prominence Configuration 

In the case of the IP configuration (Case 4b), the emerging flux rope undergoes a two-stage 
evolution as we can see from Figure 2: (i) a slow evolution (from hrs to 7 hrs) sets in upon the 
flux rope emergence; (ii) the flux rope then goes unstable at around 7 hrs (exhibiting a gradual 
acceleration) and later propagates upward into the interplanetary space with v = 70.4 km s"^. 
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The evolution of the magnetic field and plasma velocity, electric current density, and density 
enhancement for Case 4b is shown in the right columns of Figures 3, 4, and 5, respectively. It 
is clearly noted that there is no reconnection on the leading edge of the flux rope because the 
magnetic field in the upper part of the rope runs in the same sense as the external field of the 
helmet dome. However, a vertical current sheet is formed, trailing the flux rope from below as 
the rope rises (e.g.. Figs. 4g and 4h). Magnetic reconnection ensues and dissipates the current 
sheet. As a result of reconnection as well as the rising X -type neutral point due to the fixed 



boundary condition at the inner boundary jWu et al. 



19951) . the opened magnetic field is later 



reclosed low in the corona, as shown in Figures 3g and 3h. This enables the helmet streamer 
structure to recover at ~ 40 hrs (not shown). Similar processes take place in Case 4a as well 
after a vertical current sheet is formed below the old flux rope (Figs. 3c and 3d). The recovery 
to equilibrium is much slower in the IP than in the NP configuration. This is partly due to the 
different CME propagation speeds and dynamic time scales for these two magnetic topologies. It 
is interesting to note that, in the early phase, the maxima of the electric current (Figs. 4e and 4f) 
and density enhancement (Figs. 5e and 5f) occur in the vertical current sheet and the lower part 
of the emerging flux rope, respectively, in contrast with those of Case 4a. 



4. Discussions 

4.1. Operating Forces 

In a MHD representation of the corona plasma, there are three forces, viz., the Lorentz force, 
pressure gradient force, and gravitational force, which fundamentally determine the dynamics of 
the flux rope system (e.g., Wu et al. 1995; Guo et al. 1996; Hu & Liu 2000). We show in Figure 
6 the spatial distributions in the equatorial plane of the normalized radial components of these 
forces for Cases 4a (left) and 4b (right) during the early stage of the simulation. By examining 
the forces, we notice the following interesting features. 

For the NP configuration (Case 4a), the pressure force is the dominant positive force to desta- 
biUze the helmet streamer in the early stage (f < 4 hrs) (Figs. 6a and 6b), but it drops significantly 
afterwards. This is true because the emerging flux rope carries substantial mass as well as upward 
momentum into the helmet dome and a large pressure force is accordingly developed; after the 
emergence (f > 4 hrs) there is no more momentum being added into the corona except for the 
inner boundary conditions on the photosphere to maintain the background solar wind. In con- 
trast, the Lorentz force in the newly formed flux rope is relatively small. This happens for the 
following reasons. On the one hand, the emerging flux rope carries a current in a direction op- 
posite to that in the external helmet-streamer field and magnetic reconnection rapidly dissipates 
the current near the leading edge of the old flux rope around the equatorial plane. As a result, 
the current density in that region is very small and the current distribution peaks on the flanks 
of the old flux rope (Figs. 4a, 4b, and 7a through 7c). On the other hand, because reconnection 
gradually removes the constraints of the hebnet-dome's closed field ahead of the new flux rope. 
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the new flux rope's material and frozen-in magnetic field can rise and expand readily (Figs. 3a 
and 3b), leading to a decreased magnetic field strength (Figs. 8a through 8c) and further lowering 
the current density as well as the pressure gradient force inside the volume of the new flux rope. 
Being the cross product of the current density and magnetic field strength, a small Lorentz force 
(Figs. 6a through 6c) is therefore produced in the new flux rope and its vicinity, in the presence 
of the decreased current and magnetic field. 

For the IP configuration (Case 4b), the pressure gradient force and gravitational force exhibit 
similar behaviors and magnitudes as in Case 4a (Fig. 6). However, the Lorentz force tells us 
a different story. First, unlike Case 4a, the emerging flux rope in this case bears a current in 
the same direction as that of the helmet streamer above and this current produces an additional 
positive Lorentz force which can be understood in terms of the attractive force between the two 
like-signed currents. Second, there is no reconnection on the leading edge of the flux rope (cf. 
Case 4a) and the constraints from the overlying helmet-dome arcades cannot be readily removed. 
Therefore, the rise and expansion of the flux rope are suppressed by the arcades. This leads to a 
pileup of plasma (e.g.. Fig. 5e) in the flux rope during the course of the emergence and results in 
an enhanced magnetic field (Figs. 8d through 8f) and current (Figs. 4e, and 7d through 7f). For 
these two reasons the Lorentz force appears much larger in this case compared with that in Case 
4a (Figs. 6d through 6f). We also note that the Lorentz force flips its direction across the flux rope 
center due to the circulating magnetic field: this force is negative (downward) in the upper half 
of the rope, tending to hold the rising material; it switches upward in the lower half, supporting 
the dense plasma (Figs. 6e and 6f). This feature is less evident in Case 4a, as the Lorentz force 
is smaller and the topology is more complex (involving two flux ropes rather than one) in that 
case. As to the temporal evolution, the flux rope experiences two stages as we note in § 3.2. 
Early during the emergence before the center of the flux rope comes across the photosphere, the 
pressure gradient force is the most dominant to lift off the rope plasma; however, the Lorentz 
force remains negative above the center and tends to suppress the upward motion (Fig. 6e). Later 
on, after the center appears in the corona, the positive Lorentz force below it joins the pressure 
force to work against the gravity as well as the confining Lorentz force in the upper half of the 
rope. As the pressure force dies away after the emergence completion, the positive Lorentz force 
takes over the dominance. The complex interplay of these forces accounts for the slow evolution 
stage (< 7 hrs) of the flux rope which appears as a hump in the speed-time profile (Fig. 2b). After 
that, the positive Lorentz force in the lower part of the flux rope remains dominant, responsible 
for the gradual acceleration of the flux rope starting at / ~ 7 hrs, and leads to its eventual eruption. 

4.2. Magnetic Topologies and Reconnection 

Now comes the question why the flux rope erupts faster in the NP than in the IP configuration. 
We illustrate by simulation that in the NP configuration case magnetic reconnection plays a 
direct role in launching the CME (LZ02) by forming the new flux rope; this process removes 
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the closed field lines of the helmet dome as well as their constraints of the downward Lorentz 
force ahead of the new flux rope, thereby allowing the new flux rope to escape readily. Whereas, 
in the case of the IP configuration the overlying magnetic arcades tend to confine the flux rope; 
lacking mechanisms (e.g., reconnection) to remove the confinement, the flux rope thus fails to 
reach a higher initial speed. To appreciate this point, let us further examine our simulation results, 
focusing on the topological differences. 

Firstly, the reader is reminded of the physical environment specified by this numerical study, 
(i) The characteristic plasma fi at the bottom of the corona is unity, which means that the magnetic 
energy and thermal energy are comparable, and the Lorentz force and pressure gradient force are 
roughly of equal importance, (ii) A CME is initiated by introducing a flux rope to emerge through 
the photosphere into the corona, (iii) Except the magnetic polarity of the emerging flux rope, 
all the conditions are the same for the two cases arranged in a pair (NP and IP configuration), 
which implies that topology-independent forces (e.g., the pressure gradient and gravitational 
force) would behave similarly in the two cases but topology-dependent forces (e.g., the Lorentz 
force) would not. 

With these points in mind we realize that, since the flux rope emergence injects a substantial 
amount of mass, magnetic flux, electric current and upward momentum into the corona, such 
a flux rope possesses a potential to disrupt the helmet streamer, where the increasing pressure 
gradient force would play an important role especially in the early stage of the emergence (< 2 
hrs), such as in Cases 4a and 4b (Fig. 6). In the meantime, the inputted mass also induces an 
extra downward gravitational force that competes with the the upward pressure gradient force. 
The consequent evolution heavily depends on how the corona responds to such an injection 
during the first a few hours. Since the pressure force and gravitational force are similar in a pair 
of cases, we shall pay more attention on their different Lorentz force. 

In Case 4a (NP configuration), for example, the flux rope emergence drives magnetic recon- 
nection on the leading edge of the flux rope. In turn, reconnection clears the constraints ahead, 
i.e., the closed field lines of the helmet dome, and thus allows the flux rop e to rise more read- 
ily. This is a break-out situation similar to that of 



Antiochos et al 



19991). Once reconnection 



removes all the closed field lines, the flux rope is left with an open-field channel ahead in the 
helmet streamer. Note that the confining Lorentz force almost vanishes in the equatorial plane 
above the center of the new flux rope (see Figs. 6a through 6c). Therefore, the only remaining 
major constraints that would prevent the flux rope from erupting now come from the gravity. 
Driven by the large pressure gradient force, the flux rope can thus overcome the gravitational 
pull and readily escape along the open-field channel with a high initial speed. 

Nevertheless, in the IP configuration case (e.g.. Case 4b), there is no reconnection on the 
leading edge of the flux rope to remove the overlying arcades and the corona responds to the flux 
rope emergence with a negative Lorentz force above the flux rope center, much larger than its 
counterpart in the corresponding NP configuration case (e.g.. Figs. 6d through 6f). This confining 
Lorentz force, together with the gravity, overcomes part of the upward momentum and tends to 
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suppress the emergence against the pressure gradient force. Incidentally, it is worth mentioning 
that although the magnitude of the negative Lorentz force is not large in comparison with the 
more dominant gravity, when the pressure force and gravity are in competition at some critical 
point, even a small additional force would easily turn the lever of force competition to one end. In 
this sense, the negative Lorentz force plays such a crucial role. This results in the hump-shaped 
portion in the speed-time profiles for the IP configuration cases (see Fig. 2b). Losing a significant 
amount of momentum in the very early stage of the emergence, the flux rope fails to be driven to 
a high initial speed even with the positive Lorentz force coming into play later on. 

The other NP and IP configuration cases show similar behaviors of the forces as in Cases 4a 
and 4b, respectively. As the radius and emergence speed of the flux rope become smaller, less 
upward momentum is injected into the corona, and thus it is more diflicult to work against the 
gravitational pull. Depending on the energetics of the emergence, the flux rope may erupt at a 
lower speed (e.g.. Cases la and 3b) or even fail to escape (e.g.. Cases lb and 2b). 

In brief, in the NP configuration cases, the fast CMEs result from the flux rope eruptions 
mainly driven by the pressure gradient force, with little magnetic confinement; in the IP con- 
figuration cases, the slow CMEs are driven by the pressure gradient force and positive Lorentz 
force in the lower half of the flux rope, subject to the significant drags from the confining Lorentz 
force above the flux rope center This explains the distinct speed-time profiles of CMEs in the 
two topologically different types of cases. 



5. Concluding Remarks 



We have presented MHD simulations in the form of the flux rope emergence to investigate the 
relationship between magnetic topologies and two-class CMEs as suggested by LZ02. In con- 
clusion, our numerical results demonstrate that: (i) A CME-like expulsion is more readily pro- 
duced in an NP than in an IP environment, under otherwise the same conditions. This agrees 
with the results from the analytical calculat ions of magnetic energy storage in the two types of 



prominences given by 



Zhang & Low! (120031) . (ii) Early in the event, a CME originating from an 



NP environment tends to have a higher initial speed low in the corona while a CME from an 
IP environment tends to experience a slow evolution and then erupt with a lower initial speed 
and gradual acceleration, (iii) One of the ratios of the average CME speeds for t hese two types 
of mag neti c topologies is abou t 2.3, consistent with the observations reported by 



1 19761) and 



St. Cvr et al 



Gosling et al. 



jl999l) . (iv) In an NP environment, magnetic reconnection occurs on the 
leading edge of the emerging flux rope. This reconnection removes the magnetic confining force 
produced by the closed external field ahe ad of the flux rope and launches a fast CME in a manner 



similar to the Magnetic Break-out Model llAntiochos et alJl999l) . However, in an IP environment 



with reconnection absent on the leading edge and subject to the magnetic confinement from the 
overlying arcades, the emerging flux rope either fails to erupt (e.j 
in a slow CME (e.g.. Case 4b), similar to the early works given by 



Wu etal. 


1995 


) and 


Wu etal. 
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1 1997h . Reconnection, taking place in the vertical current sheet trailing the rising flux rope, is a 
passive eff'ect of the CME expulsion^. In this sense, magnetic reconnection plays a principal role 



in gen erating a fast CME in the NP configuration but does not in the IP environment dZhang et al 



20021 



On the whole, the present study qualitatively agrees well with the LZ02 scenario. However, 
for an NP configuration, the slingshot topology predicted by LZ02 is not present in our simu- 
lations and the new flux rope formed by reconnection, as clearly shown in our results, was not 
considered in their theory. The reason to account for this difi'erence may be explained as follows. 
To produce a slingshot topology (see Fig. 2 in LZ02), one must expect that reconnection occurs 
at a single X-type neutral point on top of the flux rope and as a result the pre-event, closed exter- 
nal field lines of the helmet dome reconnect with the internal field lines of the flux rope to form 
a "slingshot". In our time-dependent simulations, as soon as the emerging flux rope touches the 
overlying helmet-dome field, a current sheet forms in between and the tearing mode instability 
may set in due to the prescribed finite resistivity. This may produce many small magnetic islands 
which may rapidly coalesce to form larger islands and eventually the new flux rope. This new 
flux rope sits between the emerging rope and the helmet dome, and on the flanks are the two 
X-type neutral points where subsequent reconnection occurs. We propose that, depending on the 
reconnection development which is diflicult to predict in advance, both the sUngshot topology (as 
in LZ02) and the topology with a new flux rope formed (as in our simulations) might be the case 
in a realistic solar environment, and both topologies would be in favor of producing fast CMEs. 

It should be pointed out that, although the CMEs are generated by emerging flux ropes into 
the corona somewhat artificially in this numerical experiment, this paper is not aimed to investi- 
gate the initiation mechanisms of CMEs which have been addressed elsewhere in the literature 
(e.g., Forbes 2000; Wu et al. 2000; and references therein). However, regardless of CME initia- 
tions, the topology-dependent behavior of the CME expulsion indeed reveals the logical connec- 
tion between the magnetic topologies and the two types of CME speed-height profiles, which is 
the main goal of this study. 

The CME speeds in our simulations are systematically lower than observed values. One of 
the reasons is that we traced the center of the erupting flux rope to obtain the corresponding 
CME speed, which is an underestimate because an observed CME speed is usually measured at 
the bright CME front that gains an additional speed due to the self-expansion of the flux rope 
relative to its center. Another reason is that the initial energy content in the emerging flux rope 
is not large enough to launch a fast CME up to 800 km s This deficiency could be remedied 
by seeking low-/? MHD solutions or by upgrading this 2-D model to a 2.5-D or a full 3-D one, 
consequently increasing the energy content. These a ttempts would be beyond the scope of the 



present paper. It is worth noting that lWu et al 



1 2004ll . alternatively, have adopted a 2.5-D model 



^ Note that reconnection here does remove a small portion of the flux from the overlying, closed field in 
the helmet dome by converting it to the flux in the outer layers of the flux rope. 
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and simulated both fast and slow CMEs (at speeds comparable to observed values) in an IP 
configuration, by invoking one or a combination of the three driving mechanisms: magnetic flux 
injection, mass drainage, and additional heating. This suggests that the two-class CME speed- 
height profiles could result from a variety of mechanisms among which the scenario of LZ02 
was an initiative example. W e would expect th at even faster CMEs will be produced if an NP 



configuration is considered by 



IWu etal. 



(120041) . 

Observations, as always, will be needed to assess a theory or a numerical study. Several as- 
pects in the LZ02 theory and in our simulations are observationally testable. One prediction of the 
LZ02 theory is that magnetic reconnection, taking place either above or blow the flux rope, plays 
distinct roles in the different field topologies of the two prominence environments. Solar flares, 
which are observational manifestations of reconnection, are hence expected to exhibit distinct 
timing and morphological beha viors in asso c iation with fast and slow CMEs. S uch behaviors 



have recently been reported by 



Zhang et al. 



i 2002h and 



Zhang & Golubl (l2003l) using TRACE 



UV/EUV observations. Hard X-ray (HXR) emissions produced by accelerated particles during 
flares provide an other tool to shed ligh t on the magnetic topology of the reconnection site and 



the current sheet JSui & Holmanl2003l) . One can test whether there exist different mor phologies 



of flar es associated with the two types of CMEs, using HXR data obtained by RHESSI llLin et al 



rpt] 



2002h or the Hard X-ray Telescope onboard the Yohkoh satellite. In particular, for limb flares, 
HXR images of the loop-top (LT) source, which is presumably located near/at the reconnection 
site, can be compared with //„ images to check whether the LT occurs above or below the erupt- 
ing prominence (if any), and to see whether such an occurrence is associated with a fast or slow 
CME, respectively. 

In the case of an NP environment, a new flux rope is formed by reconnection ahead of the 
old flux rope. If the remainder of the magnetic flux in the old flux rope is significant, the CME 
would contain two flux ropes with opposite chiralities of the magnetic field. This aspect of our 
simulations may be tested by using in situ observations of the magnetic configurations in the 
interplanetary counterparts of Earth-directed CMEs, provided that the chiralities of the flux ropes 
are conserved during their interplanetary propagation. The January 10-11, 1997 magnetic cloud 
(MC) observed by the WIND spacecraft contained ^He^^/H^ abundance similar to that of the 
streamer belt material, suggesting an association between the MC and a helmet streamer In 
addition, a very cold region of exceptionally high density was detected at the rear of the MC, 
and this dense region had an unusual composition, indicating an association with the prominence 
material. This event was interpreted to be associated with a CME on January 06 , 1997 with an 
estimated speed of 450 km s ' observed by 50//0/LASCO ( Burlaga etal.ll998 1. The very cold 
region also contained a magnetic configuration, very likely of a flux rope, though its size was 
much less than the MC. We suggest that, by respectively fitting the observed magnetic field of 
the MC and the very cold region with a flux rope model, their chiralities can be determined and 
thus will provide obs ervational evidence to test our simulation results. Recently, multiple MCs 



have been reported bv lWang et al 



1 20031) . The March 03 - 05, 2001 MCs, one of the three events 
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in their study, consisted of two flux ropes with opposite chiralities, one right-handed and the 
other left-handed. They interpreted these two MCs as a consequence of two successive CMEs 
with projected speeds of 313 km s ' and 631 km s ' , respectively. However, identifying the solar 
origin of a MC has its uncertainty and it was possible that the two MCs might result from a single 
fast CME. If this is the case, such events would provide observational supports to our model. It 
will be interesting to see, statistically, whether two successive MCs are associated with a single 
CME or successive CMEs and we look forward to such kind of observational tests to check on 
our simulations. 

Ultimately, observations of the magnetic field will provide more incisive information to verify 
or reject the LZ02 theory. A study using line-of-sight MDI magnetograms is in progress in order 
to determine if there are diff'erences in the magnetic environments of the source regions of the two 
types of CMEs (M. Zhang & J. Burkepile, private communication). In the meantime, advance in 



measuring the magnetic field in the corona by polarimetric metho ds ( iLin et alJl998 



IT, in etalJ2nnnt 



ig tne rnagnetic neia in 
J2OOO : Truiillo Bueno[2001 



Jl998 : .TudgJl998 



Trujillo Bueno et alJ2002» will be able to put the theory to 



a direct test in the near future. 
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(o) Mognetic & velocity field [t=0 hrs] 
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(b) Current (nnln = -0.32, mox = 6.22) [4=0 hrs] 




Heliocentric Distonce (R,) 

Fig. 1. The initial state of the solar corona with a helmet streamer over the equator: (a) the mag- 
netic field Unes and solar wind velocity; (b) spatial distribution of the azimuthal current den- 
sity (in contours). The minimum and maximum are shown on top of Panel b in units of 
Jq = 2.29 X 10"^ A m"^ and 18 contour levels are uniformly set between these extrema. Note the 
vertical current sheet extending upward from the helmet streamer cusp. 




Fig. 2. (a) Height-time and (b) speed-time profiles of the center of the "erupting flux rope"' for 
the eight cases listed in Table 2. The solid (dotted) lines correspond to the cases with the NP (IP) 
configuration. Note that for the speed profiles in Panel b: the initial humps at f ~ 3 hrs are results 
of the flux rope emergence; the final flat portions (except for Cases lb and 2b) are extrapolations 
of the speeds evaluated at the outer boundary (r = 7.14/?^); also the final portions of curves lb 
and 2b overlap each other. 
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(a) t=5 hrs 
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(b) t=6 hrs (f) t=^5 hrs 
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(c) t=^5 hrs 
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(g) t=20 hrs 
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Fig. 3. Temporal evolution of the magnetic field and plasma flow velocity for Case 4a (NP con- 
figuration, left column) and Case 4b (IP configuration, right column). 
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(q) i=3 hrs (min = -36.5, max = 30.1) (e) t=5 hrs (mln = -0.41, max = 8.72) 
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(b) t=6 hrs (m;n = -10.9, mox=11.7) 
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(f) i=15 hrs (mln = -3.12, max=4.05) 
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(c) t=^5 hrs (nnin = -5.52, mQx = 6.62) 
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(g) t=20 hrs (min = -0.62, mQx = 6.65) 
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(d) t=25 hrs (min = -5.57, mox = 5.46) 
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(h) t=25 hrs (m;n = -0.27, mox = 6.70) 
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Fig. 4. The corresponding contours of the azimuthal electric current density (in units of Jq = 
2.29 X 10"^ A m"^) of Figure 3. The minima and maxima are shown on top of each panel and the 
intervals [min, max] are evenly divided into 38 and 34 contour levels for Cases 4a (left) and 4b 
(right), respectively. 
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(q) i=3 hrs (min = -0.04, max = 7.88) (e) t=5 hrs (min = -0.09, max=10.7) 
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(b) t=6 hrs (min = -0.32, max=4.49) 
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(f) t=^5 hrs (mIn = -0.34, max=8.68) 
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(c) t=^5 hrs (nnin = -0.30, mQx = 0.97) 
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(g) t=20 hrs (min = -0.28, mQx = 6.99) 
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(d) t=25 hrs (min = -0.28, mox = 0.73) 
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(h) t=25 hrs (min = -0.36, mox=1.49) 
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Fig. 5. The corresponding relative density enhancement, i.e., [p(r, 6, t) - p(r, 0, 0)]/p(r, 0, 0), of 
Figure 3, plotted in the same manner as Figure 4 except that here are 15 contour levels filled with 
gray scale shades (bright: low, dark: high). 
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Fig. 6. Radial components of the forces (i.e., Lorentz force , pressure gradient force , 

and gravitational force ) per unit volume as functions of the heUocentric distance in the 

equatorial plane at various times for the NP (Case 4a, left column) and IP (Case 4b, right column) 
configuration. Distance is logarithmically plotted to show details near the solar surface. All the 
forces are normalized by the gravitational force at the equator on the photosphere in the initial 
state. Note the flux rope center (i.e., the O-type neutral point) in each case. Diamond symbols 
denote the radial positions of the center of the originaUy emerging flux rope; triangles mark those 
of the new flux rope formed by recormection in the case of the NP configuration. 
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Fig. 7. Azimuthal current density distribution, (in units of 7o = 2.29 X 10"^ A m"^) at various 
times corresponding to Figure 6, in the equatorial plane {Q = 90°, solid lines) and along 6 = 72° 
(or 6 = 108°, dashed Unes) where the split current sheets in Case 4a are roughly located. The 
dotted line in each panel marks the zero level. 
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Fig. 8. 0-component of the magnetic field (Be) vs. heliocentric distance in the equatorial plane 
(0 = 90°) at various times, corresponding to Figure 6. The dotted lines mark the zero levels. Note 
that Br vanishes at = 90° under the symmetric boundary condition. 
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Table 1. Characteristic parameters of the itiitial state. 
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no (electron number density) 3.2 x 10* cm"^ 

Jo (temperature) 1.8 x 10* K 

Bo (magnetic induction strength) 2.0 G 

ySo (plasma yS) 1.0 

Cj (sound speed) 176.7 km s~' 

Va (Alfven speed) 243.9 km s"' 

Vsw (solar wind speed) 209.7 km s~' 

Note. — AU the quantities refer to the condition in the equatorial plane. The solar wind speed is evaluated 

at the outer boundary and the others at the inner boundary. 
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Table 2. Characteristics of studied cases. 



Cases 




Field<*) 


Energy*'^' 




J 






(Rs) 




(lO'i ergs) 


(kms-') 


(kms-i) 


(kms-i) 


la 


0.10 


NP 


1.76 


9.7 


174.6 


90.0 


2a 


0.15 


NP 


3.95 


14.5 


215.8 


133.9 


3a 


0.20 


NP 


7.02 


19.3 


253.6 


152.6 


4n 
Hd 




MP 


1 1 n 

L LA) 






lOl.U 


lb 


0.10 


IP 


1.76 


9.7 


0.1 


1.0 


2b 


0.15 


IP 


3.95 


14.5 


0.2 


1.5 


3b 


0.20 


IP 


7.02 


19.3 


88.0 


9.3 


4b 


0.25 


IP 


11.0 


24.2 


227.6 


70.4 



The radius of the emerging flux rope. 

Magnetic field topologies in terms of NP or IP configurations. 

The combination of the thermal and magnetic energy in the volume of the emerging flux rope, assuming 

the third-dimensional thickness Az' = O.l/Jj. 

The flux rope emergence speed in Step 2 as described in § 2. 

'•"^ The final speed of the center of the "erupting flux rope"' . This speed is evaluated at the end of the simu- 
lation (40 hrs) or at the time when the flux rope center reaches the outer boimdary of the computational 
domain, whichever occurs first. 

^ The average speed of the "erupting flux rope"' center over the time interval during which the center 
remains in the computational domain. 



